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ABSTRACT 


The  computer  codes  BLAST,  DORF,  HELP  and  HEMP  have  been  used  to 
simulate  the  flow  field  at  the  muzzle  of  a  M-16  rifle,  without  any  deflec¬ 
tors,  compensators,  or  flash  hiders.  The  flow  was  modeled  as  an  under¬ 
expanded  jet,  of  a  compressible,  inviscid  and  non-reactive  gas.  The 
calculations,  restricted  to  the  initial  stages  of  formation  of  the  jet 
were  done  on  a  cylindrically  symmetric  geometry  without  the  presence  of  a 
flow  obstruction  in  the  form  of  a  projectile.  Results  are  presented  as 
pressure  contours  as  well  as  pressure  and  velocity  profiles  at  different 
angles  from  the  centerline  of  the  muzzle  of  the  barrel,  and  plots  of  pres¬ 
sure  versus  distance  from  the  muzzle.  The  calculations  have  been  compared 
with  available  photographic  records  of  the  process. 
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I.  INTRODUCTION 


Improved  accuracy  of  fire  is  the  overriding  concern  of  evolving 
weapons  technology.  In  this  quest,  the  designer  is  greatly  aided  by  the 
availability  of  combustion  gases  which  can  compensate  for  launch  errors. 
Muzzle  devices,  by  properly  channeling  these  gases,  accomplish  this  task. 

An  extensive  survey  of  the  literature^  however,  revealed  that  no 
theoretical  analysis  exists  which  could  serve  as  a  basis  for  the  design  of 
these  devices.  Indeed,  no  gas  flow  model,  able  to  predict  the  flow 
parameters  with  any  degree  of  accuracy,  has  appeared. 

Recent  attempts  at  predicting  the  flow  field  at  the  muzzle  of  a 
gun  have  only  met  with  partial  success.  Aside  from  the  calculation  of 
overpressure  at  some  distance  from  the  gun,  such  as  those  of  references 
2  through  5,  the  codes  employed  have  been  unable  to  give  flow  details 
needed  for  rational  design  improvements.  Indeed,  some  of  the  results 
do  not  compare  favorably  with  experimental  results.  Examples  are  the 
wrong  flow  shape  in  reference  6,  the  wrong  Mach  disk  location  in  ref.  7, 
or  even  the  failure  of  the  widely  used  RIPPLE  code  to  work  beyond  a  few 
microseconds  of  real  time.® 

In  the  light  of  this  assessment  it  was  felt  that,  as  a  first  step 
toward  establishing  our  muzzle  blast  calculational  capability  and  before 
embarking  on  a  new  code  development,  the  most  promising  codes  available 
at  BRL  should  be  tried  to  ascertain  whether  they  could  ’'simulate"  the 
flow  in  fine  enough  detail  at  the  muzzle  of  a  M-16  rifle.  The  codes 
chosen  were  BLAST®,  DORF^,  HELP  ,  and  HEMP^t  This  report  summarizes 
the  undertaking. 

The  report  is  organized  as  follows:  in  chapter  two  the  flow  model 
and  the  initial  conditions  are  discussed.  This  is  followed  by  a  descrip¬ 
tion  of  the  salient  features  of  the  computer  codes  in  chapter  three.  The 
calculational  results,  including  the  criteria  for  evaluation  is  given  next. 
We  summarize  our  conclusion  in  the  last  chapter. 

II.  DESCRIPTION  OF  THE  FLOW  TO  BE  SIMULATED 


A.  The  Flow  Model 


The  dynamics  of  the  development  of  the  flow  is  illustrated  in  Figure 
lo  In  particular,  as  the  projectile  starts  to  move  down  the  gun  tube,  it 
pushes  a  column  of  air  ahead  of  it.  This  compressed  air  exits  the  tube 
pushing  ahead  of  it  a  spherical,  weak  precursor  shock.  In  the  mean  time, 
while  the  projectile  is  still  in  the  tube,  the  hot,  high  pressure  gases 
leak  around  it  and  escape  from  the  tube. 
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As  the  projectile  emerges  from  the  tube  it  has  to  traverse  and 
interact  with  the  cloud  of  burnt  gases.  At  this  stage  the  propellant 
gases  actually  move  faster  than  the  projectile.  Consequently,  in  addi¬ 
tion  to  a  precursor  shock  ahead  of  the  projectile  and  the  Mach  disk 
formed  by  the  exhaust  gases  at  a  few  calibers  from  the  tube,  we  also 
have  a  backward  facing  shock  at  the  base  of  it.  At  a  distance  of  ten 
to  fifteen  calibers  from  the  gun  tube,  the  picture  changes,  the  bullet 
having  acquired  sufficient  kinetic  energy  to  enable  it  to  outrun  the 
propellant  gas ,  punches  through  the  precursor  shock  and  enters  the  free 
flight  stage.  At  this  point  our  interest  in  the  gas  flow  ceases. 

In  the  presence  of  muzzle  attachments,  such  as  compensators  or 
flash  hiders,  the  flow  is  much  more  complicated  with  the  flow  interacting 
with  the  flew  obstruction  resulting  in  the  reflection  of  shock  waves 
and  rarefactions.  The  present  study  doesn't  consider  the  effects  of 
such  devices.  Rather  we  model  the  development  of  the  flow  at  and  near 
the  muzzle  at  early  times,  that  is,  within  the  first  few  microseconds 
of  the  process. 

The  pressure  inside  the  tube  is  much  higher  than  the  ambient  pres¬ 
sure.  As  the  gases  exhaust  through  the  muzzle,  the  jet  expands  to  the 
ambient  pressure  increasing  its  cross  sectional  area  and  creating  the 
"shock  bottle"  configuration  of  an  underexpanded  jet.  This  comes  about 
as  follows:  as  the  expansion  waves  from  the  lip  of  the  nozzle  reach 
the  constant  pressure  boundary  of  the  jet,  they  are  reflected  as  compression 
waves  and  coalesce  to  form  a  shock.  For  slightly  underexpanded  nozzles, 
the  intercepting  shocks  meet  at  the  flow  axis  leading  to  the  diamond 
configuration.  For  great  pressure  ratios,  however,  the  shocks  are  con¬ 
nected  by  a  Mach  disk.  Past  the  Mach  disk  the  flow  is  subsonic  bounded 
above  by  the  slipstream  originating  from  the  triple  point. 

B.  The  Equations  of  Motion  and  Boundary  Conditions 

The  flow  simulation  was  done  on  an  axisymmetric  geometry  as  shown  in 
figure  2.  We  assume  that  there  are  no  muzzle  attachments  and  no  flow 
obstruction,  in  the  form  of  a  projectile,  is  present.  The  state  of  the 
high  pressure  gases  was  taken  to  be  equivalent  to  that  of  the  combustion 
gases  generated  in  the  firing  of  a  M-16  Tifle  with  a  1.795  x  10“'5kg 
charge.  Since  we  were  only  trying  to  simulate  the  first  few  micro¬ 
seconds  of  the  history  of  the  flow,  the  tube  containing  the  high 
pressure  gases  was  allowed  to  empty  and  no  attempts  were  made  to  impose 
either  choking  or  real  time  dependent  conditions  at  the  muzzle.  The 
surrounding  were  taken  to  be  standard  atmosphere  at  rest  while  the 
combustion  products  were  taken  to  be  a  gas,  without  particulate  matter,  with  a 
y-1.24  at  the  conditions  given  later  in  this  section. 
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Figure  i  Stages  of  Muzzle  Blast  Development. 

a  The  Jet  Formed  by  the  Blowby  Gases  wifh  the  Projectile  still  in  the 
Gun  Barrel.  Note  the  Mach  Disc  and  the  Triple  Shocks  Surrounded 
by  a  Region  of  Intense  Turbulence. 


b.  Generation  of  the  Secondary  Blast  Due  to  the  Explosion  of  the 
Unburnt  Propellant  Gases.  An  Interesting  Feature  of  the  Flow  Field 
is  the  "Bow  Shock  "  on  the  Base  of  the  Projectile. 


Shadowgraphs  Courtesy  of  Mr.  D.  Shear,  EBL  of  BRL,  M-16  Projectile 
from  Mann  Barrel.  Velocity  of  Projectile  970  m/s. 


The  Noble-Abel's  equation  of  state  gives  a  good  approximation 
to  the  product  gases  in  the  gun  tube  and  at  the  muzzle.  It  can  be  written 

aS  RT  1  m 

P  =  1 -  FT  (1) 

-  -  6  i 
P 

where  0  is  the  co-volume,  whose  value  in  our  case  is  approximately 

1.0cm^/g>  Mi  is  the  mole  mass  and  R  the  universal  gas  constant. 

In  a  cylindrical  co-ordinate  system,  stationary  with  respect  to  the 
gun  tube,  the  conservation  laws  for  an  inviscid,  compressed  gas  are  the 
following: 

ifi.  .  3Qpu)  3(pv)  =  n  m 

3t  3r  3z  "  J 

p  .  u  2“  ♦  v  .  |E  .  0  (3) 

^  3t  3r  9z  dr  v  1 


21  .  u  2l  +  v  2l  .  iE  . 


pH*ulf'tvlT  +  I  2E™.  .  |H.  .  o  C5) 

3t  3r  3z  r  3r  3Z 

These  are  the  conservation  of  mass,  momentum  in  two  directions,  and 
energy,  respectively.  In  addition,  an  equation  of  state  must  be  specified. 

For  boundary  conditions,  the  impermeability  of  the  walls  as  well 
as  reflection  at  the  centerline  and  free  outflow  at  the  ends  or  a  rigid 
wall  at  one  of  the  boundaries  were  specified. 


The  initial  conditions  used  in  the  calculations  were  as  follows: 

1.  Gun  Tube  (Mann  barrel) 

caliber  5.56  x  10^  [®] 

thickness  of  barrel  0.50  x  10_i  [m] 

length  of  barrel  6.93  x  10" *  [m] 


2.  Propellant 

total  mass  of  propellant 
propellant  density 
flame  temperature 
explosive  force  (zero  deterrent) 
mole  mass 

co- volume 


1.795  x  10“3[kg] 
1.605  x  103 
3.360  x  103  [  °K] 
1.115  x  10s  [m] 
2-56  * 10'2 

_  3 

1.174  x  10"3[^] 
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Ambient  Conditions 
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The  calculational  grid  chosen  was  from  the  muzzle  to  a  point 
approximately  four  calibers  downstream  and  from  the  centerline  two 
calibers  radially.  In  BLAST  twenty  mesh  points  were  taken  to  be  inside 
the  barrel.  In  DORF  and  HELP,  ten  by  ten  and  in  HEMP  five  mesh  points 
were  alotted  for  the  high  pressure  gases  and  the  thickness  of  the  barrel 
was  taken  as  twenty  by  ten  except  in  BLAST  where  it  was  one  mesh  width 
thick.  The  meshes,  rectangular  in  geometry,  each  measured  0.25  x  10“^m 
on  a  side.  The  different  cell  numbers  considered  were  dictated  by 
storage  limitation  of  the  BRL  computing  facility. 
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Figure  2.  Geometry  of  the  Region  of  Interest. 


TABLE  I 


SUMMARY  OF  CODES  CONSIDERED 


Code  Name 

BLAST 

DORF 

HELP 

HEMP 

Reference 

6 

8 

9 

10 

Geometry 

Axi symmetric 

Axi symmetric 

Axi  syrrane  trie 

Axisymmetric 

Planar 

Co-ordinates 

Eulerian 

Eulerian 

Eulerian 

Lagrangian 

Number  of 
Materials 

1 

2 

15 

Multi 

Equation  of 
State 

Ideal  gas 

Ti llotson 

Ideal  gas 

Tillotson 

Wide  choice 
avai  lab  le 

Artificial 

Viscosity 

No 

Avai lab le 

Yes 

Grid 

Size 

110  x  40 

80  x  40 

40  x  40 

20  x  20 

Special 
Feature (s) 


Uses  Godunov's 
Method 


Strength  Option 


Tracer  Particles 
Strength  Option 


Requires  Large  Storage 
Strength  Option 


III.  DESCRIPTION  OF  THE  COMPUTER  CODES 


This  chapter  outlines  the  main  features  of  the  four  computer  codes 
used  in  the  simulation  of  muzzle  blast.  For  further  details  the  reader 
is  referred  to  the  cited  references. 

A.  BLAST6 


The  computer  code  BLAST,  designed  to  solve  the  inviscid,  axisymmetric 
compressible  flow  equations  of  a  one  component  gas  with  a  moving  boundary 
representing  the  projectile,  is  based  on  the  Godunov  difference  technique. 
The  advantage  of  this  method  is  that  it  can  approximate  discontinuities  with 
a  monotonic  profile,  but  on  the  other  hand,  it  is  only  first  order  accurate. 
Godunov's  scheme  consists  of  two  steps.  First  intermediate  values  are 
calculated  by  solving  the  Riemann  problem  for  the  resolution  of  a  discon¬ 
tinuity,  and  then  the  updated  values  are  obtained  by  leapfrog  differencing. 
In  one  dimension,  then  typically. 


wn+1  =  wn  _  AI  _  fH+1/2 

m  n  h  ^  m+1/2  rn-1/2"^ 


(6) 


where  w  and  f  are  vectors  and  n+1/2  represents  the  intermediate  time  level  , 
The  stability  restriction, 

O  =  r^-  (U+C)  <1  (.7) 

h  max  “ 

is  the  usual  Courant  condition. 


In  BLAST  there  is  an  innovation  to  this  technique.  Using  the  idea 
of  splitting  due  to  Yanenko,  instead  of  solving  a  two  dimensional  problem, 
two  one  dimensional  problems  are  solved.  That  is,  the  idea  is  to  integrate 
one  set  of  equations,  starting  with  a  prescribed  set  of  initial  conditions 
for  a  half  time  step  along  one  of  the  co-ordinate  directions.  Then,  using 
these  results  as  initial  conditions,  the  second  set  of  equations  is  integrated 
in  the  other  co-ordinate  direction  during  the  second  half  time  step.  One 
of  the  advantages  derived  from  this  method  is  that  proven,  high  order  explicit 
numerical  schemes  developed  for  one  dimensional  unsteady  flow  equations  can 
be  used  for  the  integration. 


To  start  the  calculation  conditions  within  the  barrel  are  prescribed 
to  be  the  same  as  those  existing  just  as  the  projectile  exits  the  barrel. 

The  gas  was  taken  to  be  ideal  and  with  y  constant.  At  the  solid  boundaries, 
the  normal  velocity  is  set  equal  to  zero.  The  downstream  boundary  is 
expanding  just  ahead  of  the  blast  wave  up  to  an  adjustable  limit.  Past 
this  point,  the  gas  is  allowed  to ’flow  through"  the  boundary.  There  is 
no  flow  allowed  across  the  center  line  of  the  tube.  The  program  contains 
the  option  of  the  presence  of  a  projectile  moving  at  a  constant  velocity. 
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B.  DORF9 


This  code  was  developed  for  hypervelocity  impact  calculations,  but  is 
general  enough  to  permit  purely  hydrodynamic  flow  simulation.  It 
consists  of  two  parts:  TCLAM,  the  generator  code  and  DORF  where  the 
actual  calculations  are  carried  out. 

Both  of  these  codes  were  written  for  cylindrical  geometry  with  the 
left  hand  boundary  being  the  axis  of  symmetry.  The  region  of  interest 
is  subdivided  into  cells  with  particles  assigned  to  each.  The  particles 
are  proportional  to  the  densities  and  energies  such  that,  for  example, 
the  total  mass  is  determined  by  summing  all  the  particles.  This  version 
of  the  code  can  handle  two  different  materials. 


TCLAM  sets  up  the  problem.  This  involves  specifying  the  geometry, 
grid  dimensions  and  initialization  of  the  variables.  In  addition,  tracer 
particles  are  defined,  if  used,  and  constants  for  the  density,  internal 
energy  and  velocity  curve  fits  are  given. 


In  DORF,  the  maximum  grid  can  not  exceed  (IMAX)  (JMAX)  <_  4200, 

Inside  this  grid  there  is  an  active  grid  consisting  of  that  part  of  the 
problem  where  there  is  actual  movement  of  mass  or  internal  energy.  The 
calculation  is  done  in  two  steps.  The  first  step  solves  the  momentum  and 
energy  equation  as  a  function  of  the  pressure  forces.  In  the  next  phase 
the  transport  terms,  omitted  from  the  first  phase  are  approximated  by 
transporting  mass,  momentum  and  energy  across  the  cell  boundaries. 


The  finite  difference  form  of  the  conservation  laws  are  as  follows: 
The  radial  momentum  equation. 


^k  pkL  ~  pkR 
k  9t  ~  2Ar 

k 


(8) 


where  the  subscript  kL  (kR)  indicates  the  state  to  the  left  (right)  of  the 
cell  under  investigation  and  Ar,  is  the  cell  size  of  the  k-th  cell  in  the 
r  direction.  Analogously,  for  the  axial  momentum 
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where 


(11) 


I 

and  n  stands  for  the  changes  due  to  pressure  forces  only.  The  left  hand 
sides  of  the  equations  are  differenced  in  the  conventional  manner,  i.e. , 


t 


.n 


Tn  „ 

Ik  etc. 


(12) 


The  solution  of  the  energy  equation  requires  the  velocities  at  two  dif¬ 
ferent  time  steps  so  that  two  passes  through  the  first  phase  of  the  cal¬ 
culation  are  required.  The  first  pass  integrates  the  momentum  equations  and 
and  using  the  old  velocities  obtains  preliminary  interface  velocities. 

The  second  phase  bypasses  the  momentum  equations  and  computes  the  new 
interface  velocities.  Details  are  given  in  reference  9. 


The  DORF  code  admits  two  types  of  boundary  conditions:  reflective 
or  transmittive .  The  left  boundary,  the  axis  of  symmetry;  is  always  re¬ 
flective,  no  flow  being  allowed  across  it.  The  top  is  always  transmittive 
while  the  bottom  can  be  either  reflective  or  transmittive. 


Upon  completion  of  the  first  phase  of  the  calculation,  involving 
the  momentum  and  energy  equations,  the  second  phase  is  executed.  Here, 
the  transport  terms  omitted  in  phase  one  are  approximated  by  transporting 
mass  across  the  cell  boundaries.  This  entails,  of  course,  transfer  of 
momenta  and  energy. 


The  conservation  of  mass 


,  in  finite  difference  form  is 

"  2rkVk 

k  L(rkL  *  rk)ark  '(rk  *  rkL)trk 


+ 


+ 


PkBvkB 

Azk 


(13) 


Weighting  of  the  velocity  in  the  mass  flux  equation  must  be  used. 
DORF  performs  this  in  the  manner  indicated  in  figure  3. 
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Assuming  flow  from  left  to  right,  the  mass  that  crosses  surface  i  is 

that  between  i  and  ii.  Let  A  =  i-ii  =  un  t  with  un  being  the  weighted 
velocity  in  A.  Upon  expanding  the  expression  for  the  mean  velocity  u^ 

“i  -  ^  ^  *  v.i  <14> 

in  a  Taylor  series  about  i-A,  the  following  expression  is  obtained 


A_ 

At 


~°-5  »  Vi) 

tvrV  +1 


(15) 


is  used  in  the  mass  flux  equation  if  in  the  above  expression 
+  up+^)  >  0,  otherwise  is  used. 


equation 


(15) 


One  can  expect  flux  of  mass,  momentum  and  energy  across  all  four 
sides  of  a  cell.  The  new  velocities  and  internal  energy  is  the  difference 
that  was  there  before  and  after  the  movement.  Also,  after  the  fluxes  have 
been  determined,  it  must  be  ascertained  how  much  of  each  material  moved. 
The  different  approaches  are  again  given  in  the  cited  reference. 


The  direction  of  sweep  for  the  numerical  integration  is  important 
because  it  can  yield  preferential  mass  treatment  depending  on  whether 
rows  or  columns  are  swept  first.  This  problem  is  eliminated  by  a  scheme 
of  "looking  ahead".  For  details  see  reference  9. 


The  time  step  for  the  integration  is  determined,  as  in  the  PIC  method, 
being  based  on  the  particle  velocity.  Again,  it  is  the  Courant  condition 
that  must  be  satisfied.  Finally,  the  equation  of  state  used  is  the  one 
due  to  Tillotson.  See  the  discussion  of  the  HELP  code  for  particulars. 


DORF  contains  an  artificial  viscosity  subroutine,  which  may  be  called, 
as  needed.  The  code  also  has  a  Lagrangian  feature  incorporating  tracer 
particles.  Their  motion  is  determined  by  using  the  velocities,  pressure  and 
energies  obtained  in  the  calculation.  Thereis  no  influence  on  the  hydrodynamics 
from  the  tracers,  but  the  tracers  give  a  good  idea  of  the  development  of  the 
flow  field. 

C.  HELP9 


HELP  is  a  two  dimensional,  axi symmetric,  Eulerian  code  developed  for 
solving  hydrodynamic  elastic-plastic  flow  problems.  The  algorithm  was 
derived  from  the  particle  in  cell  (PIC)  method  with  the  discrete  particles 
replaced  by  a  continuum  to  allow  the  treating  of  small  compressions  without 
the  need  to  resort  to  too  many  particles.  The  version  available  to  us 
was  limited  to  only  one  material. 
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The  differential  equations  expressing  the  conservation  laws  are 
replaced  by  analogous  integral  equations  obtained  by  integrating  over 
a  cell  volume  V,  followed  by  the  conversion  of  the  volume  integral  of 
divergences  to  surface  integrals  over  the  cell  surface  yielding 
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The  difference  analogs  are  explicit  and  accurate  to  terms  of  first  order. 
The  calculation  proceeds  in  phases5 

Phase  I  Acceleration  of  the  material  in  a  cell  due  to  pres¬ 
sure  and  the  deviatoric  stress  gradients  across  the  cells  are 
calculated. 

Phase  II  Updated  velocities  are  used  to  transport  the  ap¬ 
propriate  amount  of  mass,  momenta  and  internal  energy  across 
each  cell  boundary  and  the  revised  cell  quantities  are  com¬ 
puted. 

Ph as e  III  The  deviator  stresses  acting  on  each  cell  edge  and 
the  hoop  stresses  are  determined  here.  For  the  gas  phase  this 
step  is  bypassed. 


In  the  multi-material  version  of  the  code,  each  material  is  circumscribed 
by  a  series  of  massless  tracer  particles  which  are  propagated  at  the  local 
average  velocity  for  the  continuum.  At  the  boundaries  and  material  inter¬ 
faces  diffusion  is  prevented  by  employing  Lagrangian  techniques. 


HELP  uses  the  Tillotson  equation  of  state 

b 


[a 


■]  Ip  +  Ay  +  Bp 


Ion 


+  1 


where  a  =  material  constant 
b  =  material  constant 


A  =  material  constant 
B  =  material  constant 
y  =  ry-1 

I  =  specific  internal  energy 
p  =  mass  density 


(19) 


22 


A  is  equal  to  pQ  C  where  C  is  the  speed  of  sound  and  pQ  the  zero 
pressure  density,  a  is  the  asymptotic  Gruneisen  coefficient  and  b  equals 
GQ-a  with  Gq  the  zero  pressure  Gruneisen  parameter.  B  is  so  chosen  as 
to  fit  the  experimental  data. 

The  deviatoric  s tress  increments  dS  —  are  determined  by  the  elastic 
relation 

dSi =  2Gde'i  .  (20) 

•1  J 

where  G  is  the  modulus  of  rigidity  and  de^ • *  is  the  increment  in  the 
deviatoric  strain.  ^ 

For  air,  a  y-law  equation  of  state  in  the  form 

p  =  (y-l)pE  (21) 

where  E  is  the  specific  internal  energy,  is  used. 

The  code  has  a  variable  zoning  option  as  well  as  a  rezoning  feature. 
The  latter  can  be  useful  to  eliminate  oscillations  when  a  shock  passes 
through  a  transmittive  grid  boundary.  The  boundaries  require  special 
treatment  because  quantities  are  not  defined  for  cells  outside  the  grid. 

One  may  specify  a  transmittive  boundary  at  the  top  and  right  of  the  grid 
and  either  transmittive  or  reflective  boundary  at  the  axis.  At  border 
cells  boundary  conditions  are  specified  by  assuming  imaginary  cells  out¬ 
side  the  grid.  If  the  boundary  is  transmittive,  the  imaginary  cell  has 
the  same  variable  values  as  the  border  cell  while  for  a  reflective  boundary, 
the  states, save  for  a  sign  of  the  normal  velocity  component,  are  the  same. 

Stability  is  ensured  by  using  for  At  the  minimum  value,  over  all  cells, 
of  the  smallest  cell's  dimension  divided  by  the  sum  of  the  cell's  maximum 
velocity  component  and  its  sound  speed. 

D.  HEMP11 


Of  the  four  codes  tried  in  this  simulation,  only  HEMP  is  Lagrangian. 

It  was  intended  for  dynamic  plasticity  investigations,  but  has  also  been 
used  for  flow  simulation^.  This  is  an  extremely  versatile  code,  especially 
with  regard  to  treating  solid  deformations,  and  even  includes  a  combustion 
option.  A  wide  variety  of  equations  of  states  is  available  including 
those  for  gases.  The  code  includes  an  artificial  viscosity  so  that 
developing  shocks  are  treated  automatically  in  the  course  of  computation. 

Much  of  the  innovation  in  HEMP  concerns  the  treatment  of  strength  of 
materials  effects.  These  are  of  no  interest  for  our  simulation.  Below, 
we  summarize  some  of  the  salient  features  of  the  code  relevant  to  the 
hydrodynamic  calculations. 
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The  scheme  for  mass  zoning  is  shown  in  figure  4.  Basically,  the 
material  is  divided  into  quadrilaterals  with  a  grid  j-k  that  moves 
with  the  material.  The  mass  at  time  zero  of  each  quadrilateral  is  ob¬ 
tained  as  the  product  of  the  initial  density  and  volume  swept  out  by  the 
quadrilateral  rotated  about  the  x-axis.  Thus  for  quadrilateral  1 


M, 


<!  ■ 1/3  <^>i  *  *3  *  yp  Aa  * <>r  *  y2  *  r? 

Aa  is  an  area  and  is  determined  from 
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(Aa)i  =  1/2  [x2  (y3  -  y4)  -  x3  (y4  -  y2)  +  x4(y2  -  y3) ] 

and  analogously  for  A^. 

The  conservation  of  mass  then  reads  as 

Vj  -  1/3  $)  t  Uy”  *  y?  -  yA)  A"a  ♦  (yj  *  y^  *  yj)  ^  (24) 
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The  equation  of  motion  in  the  x-axial  direction  is 
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and  analogously  for  the  y-direction»  For  explanation  of  subscript  notation 
see  Figure  4, 


The  internal  energy,  E  is  determined  from 


e"+1  =  [En  -  (Pn  +  q)  (Vn+1-Vn)  +  Az 
where  Az  is  the  change  in  distortion  energy. 
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Figure  4.  The  Layout  of  the  Computational  Grid  of  the  HEMP  Code. 
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(30) 


P"+1  =  P(En+1,  Vn+1)l 


The  artificial  viscosity  q  is  given  by 
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Boundaries  are  treated  by  reflection,  using  imaginary  cells  and  the 
stability  condition  is  that  derived  by  von  Neumann  and  Richtmyer.  The 
version  of  the  code  at  our  disposal  does  not  have  the  capability  to 
rezone  automatically,  nor  to  indicate  degenerate  cells. 
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IV.  CALCULATIONS  RESULTS 


A.  Criteria  for  Evaluation 

Ideally,  computational  results  should  be  compared  with  experimental 
observations.  Due  to  the  extreme  flow  conditions  present  at  the  muzzle  of 
a  gun,  only  a  limited  amount  of  such  data  is  available  1^,14.  Essentially 
these  are  photographic  records,  including  shadowgraphs,  at  preselected 
time  intervals  or  overpressure  measurements  downstream  of  the  muzzle. 

'Therefore,  we  compare  our  results  with  photographs  of  firings  of  an 
M-16  rifle  at  several  times  paying  particular  attention  to 

1.  the  shape  of  the  flow  field 

2„  location  of  the  shocks 

3.  velocity  of  the  shocks 

In  addition,  we  compare  the  codes  among  themselves,  for  costs,  that  is  how 
much  computer  time  does  it  take  for  a  certain  real  time  to  be  simulated  and 
how  much  storage  is  required. 


As  a  further  check,  we  compare  the  computational  results  with  the  prediction 
of  blast  wave  theory^  steady  jet  theory  .  The  latter  gives  shock, 

Mach  disk  locations  and  the  geometry  of  the  jet  such  as  L/D  ratios.  Since 
tnese  are  only  valid  for  steady  flow,  the  conclusions  must  be  interpreted 
carefully. 


Under  the  conditions  of  isentropic  f low ,  Crist  et  al^, 
following  relationship  for  the  position  of  the  Mach  disk  X^, 
of  the  pressure  ratio. 
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(32) 


There  is  ambient  pressure  and  D  represents  the  sonic  nozzle  exit 
diameter.  For  M  ,  the  Mach  number  upstream  of  the  shock,  much  greater  than 

one  and  using  the  relationship  between  the  stagnation  pressure  and  Mach 
number  as  well  as  the  relationship  between  static  pressure  behind  and  before 
the  normal  shock  and  Mach  number,  one  obtains 
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Again,  it  should  be  kept  in  mind  that  the  derivations  are  for  steady 
flow  conditions  of  an  underexpanded  jet. 

As  an  order  of  magnitude  check  on  the  location  of  the  blast  wave  as  a 
function  of  time  and  energy  deposition,  we  use  the  estimate  given  by 
Sedov^.  Under  cylindrical  symmetry,  the  shock  wave  co-ordinate  R,  measured 
from  the  point  of  explosion,  as  a  function  of  time,  is  given  by 
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(54) 


where  t  is  the  elapsed  time  and  E  the  energy  per  unit  volume  deposited 
into  the  flow. 


Finally,  at  early  times  empirical  correlations 
pressure  profile  at  the  muzzle  may  be  expressed  as. 
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B.  Results 


Under  the  initial  conditions  summarized  in  Chapter  II,  the  four  codes 
were  run  up  to  five  hours  of  machine  time  or  until  the  code  got  into  trouble, 
whichever  occurred  first.  The  results  are  summarized  below. 

First,  some  general  remarks.  The  calculations  were  started  in  a 
shock  tube  fashion,  that  is,  the  end  of  a  tube  containing  high  pressure 
gases  was  exposed  to  the  ambient  at  the  time  zero.  However,  we  assumed 
the  combustion  gases  throughout  the  tube,  at  time  t=0 ,  to  possess  a 
uniform,  finite  velocity  corresponding  to  real  muzzle  conditions  at  the 
time  the  projectile  completely  clears  the  gun  tube.  This  in  contrast 
to  the  real  process  where  the  moving  projectile  in  the  tube  compresses  a 
column  of  air  ahead  of  it  leading  to  the  first  shock  observed  on  the 
photographs.  Also,  combustion  gases  leak  around  the  projectile  leading 
to  the  precursor  blast  wave. 

In  addition,  by  excluding  the  presence  of  the  projectile,  the  calcu¬ 
lated  flow  field,  of  necessity,  deviates  from  the  photographs  of  actual 
firings.  This  includes  the  rearward  facing  bow  shock  at  the  base  of  the 
projectile  and  at  early  times,  the  region  of  the  intense  shock  interaction. 

In  the  figures  the  shock  was  located  at  the  point  of  steepest  gradient 
in  the  pressure  profile.  Admittedly,  this  is  somewhat  arbitrary  in  view  of 
the  fact  (see  figure  5,  for  example)  that  the  shock  is  spread  over  up  to 
10  mesh  widths  in  the  calculations,  but  alternative  criteria,  are  no  more 
accurate. 

The  contour  plots  (figures  11-12)  suggest  the  location  of  the  starting 
vortex  at  the  lip  of  the  gun  tube.  This  is  also  apparent  from  the  velocity 
plots  (figures  8-10)  as  one  proceeds  laterally  or  radially  from  the  muzzle. 
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None  of  the  calculated  results  showed  the  presence  of  the  Mach  disk. 

In  case  of  the  BLAST  code,  for  t>l.S8ys,  due  to  the  subsonic  throat  condi¬ 
tions,  the  disk  is  not  expected  on  physical  grounds.  Also,  for  these  times, 
the  steady  state  formula(32)  would  place  it  at  a  position  to  which  the  flow 
has  not,  as  yet,  penetrated. 

L.  BI.A5T 

The  BLAST  code  was  run  out  to  400  time  cycles  corresponding  to  20.8ys 
real  time.  The  run  took  500  minutes  on  BRLESC  II, 

The  computational  results  are  illustrated  in  the  (a)  portion  of 
figures  5  to  10  where  pressure,  axial  and  radial  velocity  as  a  function  of 
distance  along  the  tube's  center  line,  at  forty  five  degrees  to  the  muzzle 
and  ninety  degrees  to  the  muzzle,  all  measured  from  the  axis  of  the  tube 
at  the  muzzle,  are  plotted.  Each  of  the  diagrams  shows  curves  corresponding 
to  the  labeled  real  times.  Figure  16  gives  the  shock  location  as  a 
function  of  time.  Curves  of  constant  pressure  at  succeeding  times  is  illustraic 

in  figure  11.  Shock  velocity  was  observed  to  vary  between  5.755  x  10^[— ] 

3  ^ 

at  3us  down  to  about  3.28  x  10  [— ]  at  8ys.  It  should  be  noted  that  the 

barrel  thickness  in  BLAST  is  only  one  mesh  width  in  contrast  to  the  other 

codes . 

As  expected,  plots  show  a  decay  in  the  flow  variable  amplitudes  as 
one  proceeds  away  from  the  muzzle.  However,  the  shape  of  the  flow  field, 
except  for  the  bow  shock  deviates  appreciably  from  the  photograph  of  the 
flow  field  given  in  figure  1.  Note  especially  the  lateral  boundary  of 
the  jet  close  to  the  tube.  Instead  of  the  expected  sidewise  expansion 
the  jet  remains  of  roughly  constant  cross  sectional  area  for  a  few  calibers 
downstream.  Also,  the  computational  results  do  not  reveal  the  presence 
of  the  Mach  disk,  triple  point  or  the  contact  surfaces.  Thus  significant 
flow  details  are  missing  from  the  computer  results.  In  view  of  storage 
and  running  time  limitations,  22  mesh  widths  were  alotted  to  the  gun 
tube  length  containing  the  high  pressure  gases.  Due  to  the  emptying  of 
the  barrel,  it  was  found  that  at  1.50|is  the  Mach  number  at  the  throat  was 
1.165  but  already  at  2.97ys  it  dropped  down  to  0.787  and  stayed  below 
1  for  all  times  thereafter.  Thus,  for  a  realistic  simulation  a  much 
larger  barrel  should  be  considered.  This,  of  course,  necessitates  the 
availability  of  a  computer  with  a  large  memory. 

The  blast  wave  location  according  to  equation  (34)  at  time  4ps 
should  have  been  at  0.60  x  10“^m  from  the  muzzle  while  BLAST  gave  its 
location  at  1.75  x  10~^  m.  Photographic  evidence  puts  it  at  0.75  x  10~"m. 

Thus  there  is  a  deviation  of  133%  between  the  calculated  and  observed 
values.  There  is^uali tative  agreement  between  these  results  and  those 
reported  by  Spong,  though  the  latter  worked  at  considerably  lower 
pressure  levels. 
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2.  DORF 


The  best  results  were  obtained  from  the  computer  code  DORF.  The 
flow  shape,  pressure  and  velocity  history,  and  flow  detail  were  the  most 
reasonable  from  these  runs.  To  simulate  8ys  real  time,  five  hours  of 
running  time  of  BRLESC  II  was  required. 

Our  results  are  illustrated  in  the  b  portion  of  figures  5  through 
11.  Again  we  plot  pressure,  axial  and  radial  velocity  as  a  function  of 
distance  along  the  gun  tube's  axis,  at  forty  five  degrees  and  at  ninety  degrees 
to  the  muzzle,  starting  at  the  muzzle.  Each  of  the  curves  corresponds  to 
the  real  time  indicated.  Figure  12  gives  pressure  contours  at  times 
corresponding  to  those  of  the  BLAST  calculations  enabling  a  comparison  to  be 
made  between  the  results  of  the  two  codes.  Note  that  due  to  the  shortness 
of  the  barrel,  as  the  slug  of  gas  clears  the  barrel  there  occurs  an 
appreciable  pressure  drop  at  the  muzzle.  This  is  noticeable  for  times 
greater  than  3.7ysec.  Finally,  figure  16  shows  the  shock  location 
as  a  function  of  time  as  calculated  by  BLAST  and  DORF. 

The  shock  velocity  was  observed  to  vary  between  3.S46  x  103[— ]  at 

3  m  ^ 

3.5ys  down, to  2.614  x  10  [—]  at  6.Sys.  Hie  other  flow  variables  showed 

analogous  attenuation.  The  contour  plots  indicate  the  formation  of  the 
starting  vortex  at  the  lip  of  the  muzzle.  A  finer  mesh  spacing  would 
be  required  to  obtain  more  detail  on  the  vortex  dynamics.  The  Mach 
number  at  the  throat  of  the  muzzle  ranged  from  Mx=0.725  at  1.59ysec 
to  Mx=0.716  at  7.36ysec. 

A  comparison  between  figures  1  and  12  shows  that  the  calculated  flow 
shape  is  quite  reasonable  and  figure  16  indicates  that  the  shock  velocity 
is  close  to  the  values  predicted  by  BLAST  but  deviates  somewhat  from 
experimental  values.  On  the  other  hand,  we  were  unable  to  show  the 
presence  of  the  Mach  disk  or  any  of  the  slip  lines  within  the  flow. 

With  only  ten  mesh  points  alloted  to  the  high  pressure  combustion 
gas  reservoir,  the  tube  emptied  too  fast  as  evidenced  in  the  contour 
plots  at  times  greater  than  3ysec.  This  accounted  for  the  Mach  number 
drop  at  the  muzzle  and  the  rapid  change  in  the  flow  variables  at  the 
muzzle. 
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3.  HELP 


Right  from  the  outset  serious  difficulties  were  encountered  in 
running  this  code.  We  considered  a  grid  of  40  x  40  points.  To  be  able 
to  run  the  version  of  the  code  available  to  us,  we  had  to  assume  that 
the  combustion  products  and  the  ambient  to  be  the  same  gas,  but  having 
different  internal  energies.  The  tube  walls  were  specified  as  composed 
of  steel  with  a  yield  strength  of 

Y  =  1.2  x  109  [^] 

°  nr 

and  rigidity  modulus  of 

G  =  8.2  x  107  [^]  . 

nr 

As  seen  from  figure  13b  the  flow  field  develops  as  expected  but 
something  strange  transpires.  Regions  of  negative  pressures  and  internal 
energy  appear  in  the  solid  and  in  the  ambient  gas.  Material  evaporates 
locally.  The  computer  also  gives  mixed  cells  at  physically  wrong  locations. 
At  certain  values  of  the  initial  internal  energy  of  the  solid,  the  gun 
tube  wall  deforms  and  can  even  ablate  and  mix  with  the  gas.  The  high 
pressure  gas  is  nowhere  near  the  state  where  such  effects  would  be  ex¬ 
pected.  Analogous  difficulties  were  encountered  at  other  installations 
with  the  code4  . 

Several  attempts  were  made  to  remedy  the  situation.  Strengths  and 
densities  several  times  that  of  steel  were  tried  but  the  results  were  the 
same.  On  examining  the  equations  and  boundary  conditions,  one  must  conclude 
that  probably  both  are  in  need  of  revision.  Particularly,  reflective 
conditions  should  be  employed  at  solid  walls,  and  the  strength  formulation 
reexamined.  That  the  shortcomings  in  the  code  have  not  been  discovered 
earlier  probably  can  be  ascribed  to  the  fact  that  HELP  has  been  developed 
for  hypervelocity  impact  problems  where  the  pressures  and  energies  are 
of  orders  of  magnitude  larger  than  here.  Since  the  runs  with  HELP  were 
terminated  after  only  .84  microseconds  real  time,  or  60  minutes  machine 
time,  no  comparison  with  experimental  data  is  possible, 

4.  HEMP 

HEMP,  the  Lagrangian  code,  was  run  both  on  a  grid  of  40  x  40  and 
20  x  20,  each  Ax  =  Ay  =  2.5  x  10~4m  the  smaller  region  requiring  64K 
storage.  The  longest  run  could  simulate  up  to  .9psec.  real  time.  The 
code  has  a  built  in  cutoff.  That  is,  if  the  determined  time  step  is 
less  than  10'4usec,  it  automatically  stops.  The  time  step  is  determined 
from  the  ratio  of  the  smallest  dimension  of  the  cell  to  the  velocity 
there  scanning  the  cells.  When  the  grid  gets  too  distorted,  as  happened 
in  our  run,  one  of  the  dimension  of  the  cells  will  usually  become  rather 
small,  leading  to  the  observed  conditions. 
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A  run  is  illustrated  in  figure  13.  The  fatal  distortion  starts  at 
the  lip  of  the  nozzle  and  gets  progressively  worse.  The  arrows  are 
proportional  to  the  velocity  vectors  and  give  an  idea  of  the  magnitudes  of 
the  flow  variables.  At  time  . 8ysec  the  shock  is  at  25  x  10"^m  from  the 
muzzle,  and  the  pressure  jump  is  approximately  10.  At  I  =  15,  J  =  4*  for 
example  near  the  lip  of  the  nozzle,  the  originally  square  cell  has  been 
distorted  into  a  narrow  rectangular  region.  At  t=.7ysec,  near  the  lower 
limit  of  the  time  step,  the  x-con^ponent  of  the  velocity  is  2196 

[— ]  while  the  y-component  is  411  [— ]  and  the  pressure  7.904  x  10^  [-U-] . 
s  s  mz 

At  the  throat  the  Mach  number,  M  =  1.38. 

Several  schemes  were  tried  to  improve  the  results  that  is,  to  delay 
the  distortion  of  the  grid  at  the  comer.  Reasoning  from  fluid  dynamics, 
where  boundary  layers  are  observed,  we  put  into  the  row  representing  the 
fluid  next  to  the  inside  wall  of  the  gun  and  extending  out  into  the 
ambient  a  viscosity  much  greater  than  that  assigned  to  the  other  gases. 

It  was  our  expectation  that  this  would  slow  the  distortion  process.  Indeed, 
this  was  also  observed,  but  it  could  only  prolong  the  run  by  a  few  micro¬ 
seconds  real  time .  The  value  of  the  artificial  viscosity  coefficient  used 

-  4  N 

was  60  while  the  maximum  value  of  q  was  1.987  x  10  [ — .  Elsewhere  in 

the  flow  field,  the  viscosity  coefficient  was  kept  at  ^,0. 

Since  the  computational  times  reached  were  less  than  one  microsecond 
before  the  code  exited  no  comparison  with  our  photographic  records  was 
possible.  Possibly  through  a  rezoning  scheme,  the  runs  could  be  extended 
in  duration,  but  it  is  our  opinion,  that  there  is  no  rational  ground  to 
believe  that  runs  of  the  order  of  tens  of  microseconds  real  time  could  be 
achieved.  To  reach  0.8  microseconds  real  time,  it  took  a  run  of  48 
minutes  on  BRLESC  II. 


V.  CONCLUSIONS 

Under  identical  initial  conditions,  the  computer  codes  BLAST,  D0RF, 
HELP,  and  HEMP  have  been  tested  for  capability  of  muzzle  blast  flow  field 
simulation  of  a  small  bore  weapon.  It  has  been  found  that  HELP  and  HEMP 
are  unsuitable  for  the  task. 

BLAST  is  able  to  simulate  the  gross  features  of  the  flow  including 
the  blast  wave.  However,  the  shape  of  the  predicted  flow  field  for  real 
times  less  than  20ys  deviates  in  shape  appreciably  from  the  observed. 

Also,  the  code  is  unable  to  show  the  Mach  disk  and  the  triple  shocks  known 
to  be  present. 


*1  and  J  are  the  grid  co-ordinates. 
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DORF  can  also  furnish  the  global  features  of  the  flow  field  without 
resolving,  for  the  mesh  size  chosen,  the  finer  features  of  the  flow.  Both 
DORF  and  BLAST  exhibit  order  of  magnitude  agreement  with  the  theoretical 
predictions  of  the  blast  wave  velocity  and  pressure  jump  based  on  our 
criteria  given  in  Chapter  IV.  Overall,  DORF  is  much  better  than  BLAST. 
However,  to  be  useful  for  our  purposes,  DORF  would  have  to  modified  to 
enable  it  to  handle  a  projectile,  that  is  a  moving  boundary,  as  well  as 
sharpen  the  shock  profiles.  Such  changes  are  major  and  would  require 
the  expenditure  of  considerable  time  and  effort.  A  new  code  development, 
possibly  along  the  lines  of  the  flux  corrected  transport  algorithm,  appears 
to  be  more  promising  and  is  the  approach  that  should  be  pursued. 

In  conclusion  then,  we  feel  that  none  of  the  examined  codes  are  good 
enough  as  a  design  aid  in  the  advancement  of  muzzle  device  technology. 
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PRESSURE  [■&] 


[xl0'2m] 

(a) 


Figure  5.  Pressure  Versus  Distance  along  the  Axis  of  the  Gun  Tube. 

a.  Results  from  the  BLAST  Code . 

b.  Results  from  the  DORF  Code  . 


PRESSURE 


AMBIENT  CONDITION: 


LINE  OF  GUN  TUBE  [  x  >0-2m] 

(a) 

Figure  6.  Pressure  Versus  Distance  at  90°  from  the  Axis  of  the  Gun 

Tube  Measured  from  the  Muzzle. 

a.  Results  from  the  BLAST  Code. 

b.  Results  from  the  DORF  Code. 
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PRESSURE 


(  a  ) 

Figure  7.  Pressure  Versus  Distance  at  45°  from  the  Axis  of  the  Gun  Tube. 

a.  Results  from  the  B6AST  Code. 

b.  Results  from  the  DORF  Code. 
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AXIAL  VELOCITY  [■?] 
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Figure  8. 
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OF  GUN  TUBE  [  *10t2id] 


Axial  Velocity  Versus  Distance  Along  the  Axis  of  the  Gun  Tube, 
a,.  Results  from  the  BLAST  Code, 
b;  Results  from  the  DORF  Code. 
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AXIAL  VELOCITY  [™] 


OF  GUN  TUBE  [xl0'2m] 

(  a  ) 


Figure  9.  Axial  Velocity  Versus  Distance  at  45°  from  the  Axis  of  the 

Gun  Tube. 

a.  Results  from  the  BLAST  Code. 

b.  Results  from  the  DORF  Code. 
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0  .5  1.0  1.5 

DISTANCE  AT  45°  FROM  CENTERLINE  OF 
GUN  TUBE  [x10‘2m] 

(a  ) 

Figure  10.  Radial  Velocity  Versus  Distance  at  45°  from  Axis  of  the 

Gun. 

a.  Results  from  the  BLAST  Code  . 

b.  Results  from  the  DORF  Code. 
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Figure  11.  Pressure  Contour  Plots  from  the  BLAST  Code. 
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RADIAL  DISTANCE 
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RADIAL  DISTANCE 
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AXIAL  DISTANCE  ALONG  THE  GUN  TUBE  CENTERLINE 
MEASURED  FROM  THE  MUZZLE 

[  x  10 m  ] 


Figure  14.  Pressure  Versus  Distance  Along  the  Axis  of  the  Gun  Tube 
Calculated  by  HELP. 
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AXIA  L  DISTANCE  [*  I0'2  m] 

Figure  15.  The  Distortion  of  the  Lagrangian  Grid  as  Calculated  by 

HEMP  Code  [6  figures). 


52 


RADIAL  DISTANCE  10 


TIME  =  0.10058  MICROSECONDS 
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RADIAL  DISTANCE  *10 


TIME  =  0.20037  MICROSECONDS 
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TIME  =  0.30013  MICROSECONDS 
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O  MEASURED  FROM  PHOTOGRAPHS 


Figure  16.  Comparison  of  Calculated  and  Observed  Shock  location  Downstream  of  the  Muscle  as  a  Function 
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EMPIRICAL  CORRELATION  (Ref.  1) 
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Figure  17.  Comparison  of  an  Empirical  Correlat 

as  a  Function  of  Time 


for  t  <  0.4  ms 


acknowledgment 


The  author  wishes  to  thank  Ms.  Viola  Woodward,  Messrs.  Robert  Karpp 
and  John  Harjrison  for  their  help  in  running  the  computer  codes.  The 
contour  plotting  routine  used  with  the  DORF  code  was  developed  by- 
Mr.  George  Hartwig  of  AML. 


60 


REFERENCES 


1.  Zoltani,  C.K.,  "Muzzle  Devices,  A  State  of  the  Art  Survey.  Vol.  II: 
Theoretical  Approaches  to  the  Calculations  of  Muzzle  Blast."  BRL 
Memorandum  Report  No.  2294,  1973. 

2.  Schling,  W.  ,  Sowers,  D. ,  Rogers,  T. ,  "Muzzle  Brake  Study  for  MFR 
Program".  TRW  Systems  Group  Report  on  Contract  DAAF03-72-C-0018, 
September  15,  1972. 

3.  Hoo,  G.S.,  Vagla,  J.  J.,  "Use  of  a  Conical  Muzzle  Device  to  Control 
Gun  Blast".  NWL  Technical  Report  TR-2793,  August  1972. 

4.  Moore,  G.  R. ,  "Finite  Difference  Calculations  of  the  Free  Air  Gun 
Blast  About  the  Muzzle  of  a  5/54  Naval  Gun".  NWL  Technical  Report 
TR-2794,  September  19  72. 

5.  Walther,  M.  F.,  "Gun  Blast  From  Naval  Guns".  NWL  Technical  Report 
TR-2733,  August  1972. 

6.  Taylor,  T.D.,  "Calculation  of  Muzzle  Blast  Flow  Fields".  Technical 
Report  4155,  Picatinny  Arsenal,  Dover,  N.J.,  (1970). 

7.  Sinha,  R. ,  Zakkay,  V.,  Erdos,  J. ,  "Flowfield  Analysis  of  Plumes  of 
Two  Dimensional  Underexpanded  Jets  by  a  Time  Dependent  Method".  AIAA 
Journal  9,  2363-2370  (1971). 

8.  Private  Communication  form  Ma j .  C.  W.  Nelson,  IBL  of  BRL. 

9„  Payton,  D.N.,  "The  DORF  Code  User's  Manual".  AFWL-TR- 70-60 ,  Kirtland 
Air  Force  Base,  N.M.,  (1970). 

10.  Hageman,  L.J. ,  Walsh,  J.  M. ,  "HELP,  A  Multi-Material  Eulerian  Pro¬ 
gram  for  Compressible  Fluid  and  Elastic-Plastic  Flows  in  Two  Space 
Dimensions  and  Time".  Systems,  Science,  and  Software  Report 
35R-350  (1970). 

p 

11.  Wilkins,  M,  L. ,  "Calculation  of  Elastic  Plastic  Flow".  UCRL-7322, 

Rev.  I,  University  of  California,  Livermore,  (1969). 

12.  Wilkins,  M.L.,  "Calculation  of  a  Slender  Body  Moving  Through  Air  at 

Supersonic  and  Subsonic  Velocities."  AIAA  Journal  9_,  2073-2075,  (1971). 

13.  Brode,  H.L.,  Enstrom,  J.E.,  "Interior  Ballistics  and  Gun  Flash  and 
Smoke".  RAND  Corp.  Memorandum  RM-6127-PR,  October  1969. 

14.  Photographs,  courtesy  of  EBL,  BRL. 

15.  Sedov,  L.I.,  "Similarity  and  Dimensional  Methods  in  Mechanics". 
Academic  Press,  Inc.,  New  York,  1959. 


61 


16.  Crist,  S.,  Sherman,  P.  M. ,  Glas,  D.  R. ,  "Study  of  the  Highly  Under¬ 
expanded  Sonic  Jet".  AIM  Journal  £,  68-71  (1966). 

17.  Spong,  E  ,D. ,  "The  Numerical  Solution  of  the  Expulsion  of  a  Shock 
Wave  from  the  Open  End  of  a  Cylindrical  Shock  Tube".  ScD  Thesis, 
Washington  University,  St.  Louis,  MO  1972. 

18.  Private  Communication  from  Dr.  A.K.R.  Celmins  of  AML,  BRL. 

19.  Private  Communication  from  Mr.  B.L.  Morris,  MERDC,  Fort  Belvoir,  VA. 


62 


DISTRIBUTION  LIST 


No.  of 

Copies  , Organization 

12  Commander 

Defense  Documentation  Center 
ATTN:  TIPCR 
Cameron  Station 
Alexandria,  Virginia  22314 

1  Director 

Defense  Nuclear  Agency 
Washington,  DC  20305 

1  Commander 

U.S.  Army  Materiel  Command 
ATTN:  AMCDL 
5001  Eisenhower  Avenue 
Alexandria,  Virginia  22304 

1  Commander 

U.S.  Army  Materiel  Command 
ATTN:  AMCDMA,  MG  J.R.  Guthrie 
5001  Eisenhower  Avenue 
Alexandria,  Virginia  22304 

1  Commander 

U.S.  Army  Materiel  Command 
ATTN:  AMCRD,  MG  S.  C.  Meyer 
5001  Eisenhower  Avenue 
Alexandria,  Virginia  22304 

1  Commander 

U.S.  Army  Materiel  Command 
ATTN:  AMCRD,  Dr .J .V. R. Kaufman 

5001  Eisenhower  Avenue 
Alexandria,  Virginia  22304 

1  Commander 

U.S.  Army  Materiel  Command 
ATTN:  AMCRD -TE 
5001  Eisenhower  Avenue 
Alexandria,  Virginia  22304 

1  Commander 

U.S.  Army  Materiel  Command 
ATTN:  AMCRD -TP 
5001  Eisenhower  Avenue 
Alexandria,  Virginia  22304 


No.  of 

Copies  Organization 

1  Commander 

U.S.  Army  Materiel  Command 
ATTN:  AMCRD-MT 
5001  Eisenhower  Avenue 
Alexandria,  Virginia  22304 

1  Commander 

U.S.  Army  Materiel  Command 
ATTN:  AMCRD -W 
5001  Eisenhower  Avenue 
Alexandria,  Virginia  22304 

1  Commander 

U.S.  Army  Aviation  Systems 
Command 

ATTN:  AMSAV-E 
12th  and  Spruce  Streets 
St.  Louis,  Missouri  63166 

1  Director 

U.S.  Army  Air  Mobility  Research 
and  Development  Laboratory 
Ames  Research  Center 
Moffett  Field,  California  94035 

2  Commander 

U.S.  Army  Electronics  Command 
ATTN:  AMSEL-RD 
AMSEL-CE 

Fort  Monmouth,  New  Jersey 
07703 

2  Commander 

U.S.  Army  Missile  Command 
ATTN:  AMSMI-R 

AMSMI-RBL 

Redstone  Arsenal,  Alabama 
35809 

5  Commander 

U.S.  Army  Missile  Command 
ATTN:  AMSMI-RDK 

Mr.  R.  Becht  (4  cys) 
Mr.  R.  Deep 

Redstone  Arsenal,  Alabama  35809 


63 


DISTRIBUTION  LIST 


No.  of 

Copies  Organization 

1  Commander 

U.S.  Army  Tank  Automotive 
Command 

ATTN:  AMSTA-RHFL 
Warren,  Michigan  48090 

2  Commander 

U.S.  Army  Mobility  Equipment 
Research  8  Development  Center 
ATTN:  Tech  Docu  Cen,  Bldg.  315 
AMSME-RZT 

Fort  Belvoir,  Virginia  22060 

1  Commander 

U.S.  Army  Armament  Command 
Rock  Island,  Illinois  61202 

5  Commander 

U.S.  Army  Frankford  Arsenal 
ATTN:  SARFA-J8000 

Mr.  J.  Mitchell 
SARFA-K3000 
Mr.  S.  Fulton 
SARFA-J6300 

Mr.  S.  Hirshman 
SARFA-J3000 

Mr.  A.  Cianciosi 
SARFA-L4100 

Mr.  C.  Fleischer,  Jr. 
Philadelphia,  Pennsylvania 
19137 

7  Commander 

U.S.  Army  Picatinny  Arsenal 
ATTN :  SARPA-AD 

Mr.  S.  Wasserman 
Mr.  Lindner 
SARPA-FR 
Mr.  A.  Loeb 
Mr,  D.  Mertz 
Mr.  E.  Walbrecht 
Mr.  S.  Verner 
SARPA-VE,  Dr.  Kaufman 
Dover,  New  Jersey  07801 


No.  of 

Copies  Organization 

1  Commander 

U.S.  Army  Watervliet  Arsenal 
Watervliet,  New  York  12189 

1  Commander 

U.S.  Army  Safeguard  Systems 
Command 

Huntsville,  Alabama  35804 
1  Director 

U.S.  Army  Advanced  Materiel 
Concepts  Agency 
2461  Eisenhower  Avenue 
Alexandria,  Virginia  22314 

1  Commander 

U.S.  Army  Harry  Diamond 
Laboratories 
ATTN:  AMXDO-TI 
Washington,  DC  20438 

1  Commander 

U.S.  Army  Materials  and 
Mechanics  Research  Center 
ATTN:  AMXMR-ATL 
Watertown,  Massachusetts 
02172 

1  Commander 

U.S.  Army  Natick  Laboratories 
ATTN:  AMXRE,  Dr.  D.  Sieling 
Natick,  Massachusetts  01762 

1  Deputy  Assistant  Secretary  of 
the  Army  (RSD) 

ATTN:  Mr.  C.  L.  Poor 
Washington,  DC  20310 

1  Commander 

U.S.  Army  Research  Officer 
ATTN:  CRD-AA-EH 
Box  CM,  Duke  Station 
Durham,  North  Carolina  27706 


64 


DISTRIBUTION  LIST 


No.  of 

Copies  Organization 

1  Director 

U.S.  Army  Advanced  Ballistics 
Missile  Defense  Agency 
P.  0.  Box  1500 
Huntsville,  Alabama  35809 

3  Commander 

U.S.  Naval  Air  Systems  Command 
ATTN:  AIR-604 
Washington,  DC  20360 

3  Commander 

U.S.  Naval  Ordnance  Systems 
Command 

ATTN:  ORD-9132 
Washington,  DC  20360 

2  Commander  and  Director 
U.S.  Naval  Ship  Research 

and  Development  Center 
ATTN:  Tech  Lib 

Aerodynamic  Lab 
Washington,  DC  20007 

3  Commander 

U.S.  Naval  Weapons  Center 
ATTN:  Code  753,  Tech  Lib 
Code  50704, 

Dr.  W.  Haseltine 
Code  3007,  Mr.  A.  Rice 
China  Lake,  California  93555 

4  Commander 

U.S.  Naval  Ordnance  Laboratory 
ATTN:  Code  031,  Dr.  K.  Lobb 
Code  312,  Mr.  R.  Regan 
Mr.  S.  Hastings 
Code  730,  Tech  Lib 
Silver  Spring,  Maryland  20910 

3  Director 

U.S.  Naval  Research  Laboratory 
ATTN:  Tech  Info  Div 

Code  7700,  Dr.  A.  Kolb 
Code  7720,  Dr.  E.  McClean 
Washington,  DC  20390 


No.  of 

Copies  Organization 

1  Commander 

U.S.  Naval  Weapons  Laboratory 
ATTN:  Code  GX,  Dr.  W.  Kemper 
Dahlgren,  Virginia  22448 

2  Commander 

U.S.  Naval  Ordnance  Station 
Indian  Head,  Maryland  20640 

2  ADTC  (ADBPS-1 2) 

Eglin  AFB 
Florida  32542 

1  AFATL  (DLR) 

Eglin  AFB 
Florida  32542 

1  AFATL  (DLRD) 

Eglin  AFB 
Florida  32542 

1  AFATL  (DLRV) 

Eglin  AFB 
Florida  32542 

2  AFATL  (DLRA,  Mr.  F.  Burgess; 
Tech  Lib) 

Eglin  AFB 
Florida  32542 

1  AFWL  (DEV)  (WLIL) 

Kirtland  AFB 
New  Mexico  .87117 

1  ARD  (ARIL) 

Wright-Patterson  AFB 
Ohio  45433 

1  ARL 

Wright-Patterson  AFB 
Ohio  45433 

1  ASD  (ASBEE) 

Wright-Patterson  AFB 
Ohio  45433 


65 


DISTRIBUTION  LIST 


No.  of  No.  of 

Copies  Organization  Copies  Organization 


1  RTD  (RDFE,  Mr.  Sedderke) 
Wright-Patterson  AFB 
Ohio  45433 

1  Director 

National  Bureau  of  Standards 
ATTN:  Tech  Lib 
U.S,  Department  of  Commerce 
Washington,  DC  20234 

1  Headquarters 

National  Aeronautics  and 
Space  Administration 
ATTN:  Code  EP,  M.  Adams 
Washington,  DC  20546 

1  Director 

NASA  Scientific  and  Technical 
Information  Facility 
ATTN:  SAK/DL 
P.  0.  Box  33 

College  Park,  Maryland  20740 

1  Director 

Jet  Propulsion  Laboratory 
ATTN :  Tech  Lib 
4800  Oak  Grove  Drive 
Pasadena,  California  91103 

2  Director 

National  Aeronautics  and 
Space  Administration 
George  C.  Marshall  Space 
Flight  Center 
ATTN:  MS- I,  Lib 

R-AERO-AE,  Mr.  A.  Felix 
Huntsville,  Alabama  35812 

1  Director 

National  Aeronautics  and 
Space  Administration 
Langley  Research  Center 
ATTN:  MS  185,  Tech  Lib 
Langley  Station 
Hampton,  Virginia  23365 


1  Director 

National  Aeronautics  and 
Space  Administration 
Lewis  Research  Center 
ATTN:  MS  60-3,  Tech  Lib 
21000  Brookpark  Road 
Cleveland,  Ohio  44135 

1  Advanced  Technology  Laboratory 
ATTN:  Dr.  J.  Erdos 

400  Jericho  Turnpike 
Jericho,  New  York  11753 

2  ARO,  Inc. 

ATTN:  Tech  Lib 
Arnold  AFS 
Tennessee  37389 

1  Technical  Director 

Colt  Firearms  Corporation 
150  Hyshore  Avenue 
Hartford,  Connecticut  14061 

1  General  Electric  Corporation 
Armaments  Division 
ATTN:  Mr.  R.  Whyte 
Lakeside  Avenue 
Burlington,  Vermont  05401 

1  Winchester-Western  Division 
Olin  Corporation 
ATTN:  Mr.  D.  Mettill 
New  Haven,  Connecticut  06504 

1  Sandia  Corporation 

ATTN:  Aerodynamics  Dept 

Org  5620,  Mr.  R.  Maydew 
P.  O.  Box  5800 

Albuquerque,  New  Mexico  87115 

1  California  Institute  of  Technology 
Aeronautics  Department 
ATTN:  Prof.  H.  Liepmann 

Pasadena,  California  91102 


66 


DISTRIBUTION  LIST 


No .  of 

Copies  Organization 

1  Guggenheim  Aeronautical 

Laboratory 

California  Institute  of 
Technology 
ATTN:  Tech  Lib 
Pasadena,  California  91104 

3  Cal  span  Corporation 
ATTN:  Mr.  J.  Martin 
Dr.  W.  Wurster 
Dr.  G.  Skinner 
P.  0.  Box  235 
Buffalo,  New  York  14221 

2  Franklin  Institute 
ATTN:  Dr.  Carfagno 

Dr.  Wachtell 
Race  £  20th  Streets 
Philadelphia,  Pennsylvania 
19103 

1  Director 

Applied  Physics  Laboratory 
The  Johns  Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring,  Maryland  20910 

1  Massachusetts  Institute  of 
Technology 

Department  of  Aeronautics 
and  Astronautics 
ATTN:  Tech  Lib 
77  Massachusetts  Avenue 
Cambridge,  Massachusetts  02139 

1  Director 

Guggenheim  Aerospace  Labs 
New  York  University 
New  York  Heights 
New  York,  New  York  10053 


No.  of 

Copies  Organi zation 

1  Ohio  State  University 

Department  of  Aeronautics  and 
Astronautical  Engineering 
ATTN:  Tech  Lib 
Columbus,  Ohio  43210 

1  Polytechnic  Institute  of 
Brooklyn 
Graduate  Center 
ATTN:  Dr.  G.  Morretti 
Farmingdale,  New  York  11735 

1  Director 

Forrestal  Research  Center 
Princeton  University 
Princeton,  New  Jersey  08540 

1  Princeton  University 
Forrestal  Laboratories 
ATTN :  Dr .  S .  I .  Cheng 
Princeton,  New  Jersey  08540 

2  University  of  Mechigan 
Department  of  Aeronautical 

Engineering 
ATTN:  Dr.  A.  Kuethe 
Dr.  M.  Sichel 
East  Engineering  Building 
Ann  Arbor,  Michigan  48104 


Aberdeen  Proving  Ground 

Ch,  Tech  Lib 
Marine  Corps  Ln  Ofc 
Cmdr,  USALWL 
Cmdr,  USATECOM 

ATTN:  AMSTE-BE,  Mr.  Morrow 
AMSTE-TA-R,  Mr.  Wise 
Dir,  USAMSAA 


67 


Intentionally  left  blank. 


Lassified 

Security  Classification 


DOCUMENT  CONTROL  DATA  -  R  &  D 

/■Security  cfaeisifii  efion  ot  title,  body  of  ahsfrac  t  and  indexing  annotation  must  be  entered  when  /he  overall  report  lo  clasalfled) 

v  originating  activity  (Corporate  author) 

2 m.  REPORT  SECURITY  CLASSIFICATION 

USA  Ballistic  Research  Laboratories 

Unclassified 

Aberdeen  Proving  Ground,  Maryland  21005 

26.  GROUP 

N/A 

3  REPORT  TITLE 

Evaluation  of  the  Computer  Codes  BLAST,  DORF,  HELP  and 

HEMP  for  Suitability  of 

Underexpanded  Jet  Flow  Calculation. 

4.  DESCftiP  TIVE  NOTES  ( Typ e  of  report  and  Inclusive  dates) 


s  AuTHORist  (First  name,  middle  Initial,  la  at  name) 

Csaba  K.  Zoltani 


0  REPORT  DATE 

AUGUST  1973 


8»  CONTRACT  OR  GRANT  NO. 

b.  project  no.  1J562604A607 


la.  TOTAL  NO-  OF  PAGES  lb.  NO.  OF  REFS 

67  19 


9a.  ORIGINATOR’S  REPORT  NUW6ERI5) 

BRL  REPORT  NO.  1659 


10.  DISTRIBUTION  STATEMENT 


»£>.  OTHER  REPORT  NO<S)  (Any  other  numbera  that  may  be  aeaioned 
thle  report) 


Approved  for  public  release;  distribution  unlimited. 


11  SUPPLEMENTARY  NOTES  12-  SPONSORING  MILITARY  AOTiViTV 

US  Army  Materiel  Command 
5001  Eisenhower  Avenue 
Alexandria,  Virginia  22304 

13  ABSTRACT 

The  computer  codes  BLAST,  DORF,  HELP  and  HEMP  have  been  used  to  simulate  the 
flow  field  at  the  muzzle  of  a  M-16  rifle,  without  any  deflectors,  compensators,  or 
flash  hiders.  The  flow  was  modeled  as  an  underexpanded  jet,  of  a  compressible, 
inviscid  and  non-reactive  gas.  The  calculations,  restricted  to  the  initial 
stages  of  formation  of  the  jet  were  done  on  a  cylindrically  symmetric  geometry 
without  the  presence  of  a  flow  obstruction  in  the  form  of  a  projectile.  Results 
are  presented  as  pressure  contours  as  well  as  pressure  and  velocity  profiles  at 
different  angles  from  the  centerline  of  the  muzzle  of  the  barrel,  and  plots 
of  pressure  versus  distance  from  the  muzzle.  The  calculations  are  compared  with 
available  photographic  records  of  the  process. 


.1473 


R«FL4C«f  OO  ROHM  1479.  1  JAN  «4.  WHICH  IS 
OBtOLtTt  FOR  ARMY  UBS.  ££ 


Unclas 


ecurity  Classification 


Security  Classification 


Muzzle  Blast 
Gas  Dynamics 
Shock  Waves 
Numerical  Analysis 
Computer  Codes 


